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Bending of worm-like polymers carries an en- 
ergy penalty which results in the appearance 
of a persistence length Ip such that the poly- 
mer is straight on length scales smaller than Ip 
and bends only on length scales larger than this 
length. Intuitively, this leads us to expect that 
the most probable value of the local curvature of 
a worm-like polymer undergoing thermal fluc- 
tuations in a solvent, is zero. We use simple 
geometric arguments and Monte Carlo simula- 
tions to show that while this expectation is in- 
deed true for polymers on surfaces (in two di- 
mensions), in three dimensions the probability 
of observing zero curvature anywhere along the 
worm-like chain, vanishes. 

1 Introduction 

Semi-flexible polymers are usually modeled as 
space curves (either continuous or discrete) that 
possess bending rigidity. Following the linear 
theory of elasticity of slender rods, the bending 
energy of such worm-like polymers is taken to 
be a quadratic functional of the local curvature 
[1]. The curvature is a local geometrical prop- 
erty of a continuous curve - the reciprocal of 
the radius of the tangential circle at any point 
along the curve. In discrete models the curva- 
ture is the angle (per unit length) between ad- 
jacent segments, in the limit when the length of 
these segments goes to zero (see below). The 
calculation of the statistical properties of such 
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polymers involves generating a representative 
statistical ensemble of different conformations 
of space curves and weighting them by appro- 
priate Boltzmann factors. Since these statisti- 
cal weights are Gaussian functions of the local 
curvature, one expects that the most probable 
value of the curvature is zero, i.e., that the prob- 
ability distribution of curvature is peaked at the 
origin. This expectation agrees with our intu- 
ition about semi-flexible polymers: since they 
deviate from straight lines only on length scales 
exceeding their persistence length, the local cur- 
vature of a typical configuration is expected to 
vanish almost everywhere. As we will show in 
the following using both analytical arguments 
and Monte Carlo simulations, the above con- 
clusions hold only in two dimensions; surpris- 
ingly, the most probable value of the curvature 
in three dimensions is not zero! 

2 Curves, curvature and 
measure 

Consider a space curve defined by a vector r(s) 
(with respect to some space-fixed coordinate 
system), where the parameter s is the arclength. 
The shape and the orientation of the curve are 
completely determined (up to unifrom transla- 
tion) by the direction of the tangent to the curve 
at every point along its contour, t{s) = dr/ds: 

r{s) = f ds'i{s') (1) 

J 

Since the tangent is a unit vector, \i{s)\ = 1, it 
can be represented as a point on a sphere of unit 
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radius; the coordinates of this point are given by 
the two spherical angles 9 and 0. Finally, one 
can introduce a local property of the curve, the 
curvature, defined as : 







di 






ds 



(2) 



Let us introduce a discrete version of the con- 
tinuous curve of length L, by dividing it into 
segments of length As (A^As = L), such that 
the orientation of the n'th segment connecting 
points s and s + As is defined by the direction 
of the tangent at point s. One can define the 
orientation of the n'th tangent vector (i.e., the 
n'th segment) in terms of the (n — l)'th one 

tn = cos(A^)t„_i + cos(A(y9) sin(A^^)n„_i 
+ sm{Aip) sin(A6')6„_i 

(3) 

As shown in Fig. [l| AO E [0..it] is the angle 
between the n'th and the (n — l)'th segments 
and Aip E [—tt..7i] is the polar angle between the 
two planes defined by t„_2 x tn-i and tn-i x t„. 
In the continuum limit Aip becomes the angle 
of rotation of the normal n and the binormal 
6, with respect to the tangent t to the curve. 
In differential geometry this triad t, n and b is 
referred to as the Frenet frame and its rotation 
as one moves along the curve is governed by the 
Frenet-Serret equations. As we will show in the 
following, in the limit of a continuous curve the 
angle A^ is proportional to the curvature (one 
can also show that the torsion is proportional 
to the angle Aip [2]). 

The calculation of statistical averages involves 
integration over all possible configurations of 
the polymer with respect to the measure 
D{f{s)} or, alternatively, over D{t{s)}. Dis- 
cretization allows us to replace this measure by 
the product of contributions of its segments, 
D{t{s)} Iln{dtn), where the measure of a 
single segment is given by 

di = sm{Ae)d{Ae)d{A^) (4) 

The local curvature of a discrete curve can be 
defined by replacing the differentials by differ- 
ences in Eq. |2j 

AsKn = \j{in+l - in? = ^2(1 - COs(A^^)) (5) 



where we used the equalities (tn+i)^ = (^n)^ = 
1 and tn+i ■ tn = cos{A9). This definition of 
the curvature is consistent with the one used in 
the Frenet-Serret theory only if kAs << 1, in 
which case kAs ^ \A6\. Using Eqs|4]and[5]it is 
straightforward to show that: 

di{s) = d{Aip){s)As^K{s)dK{s) (6) 

The above discussion applies to curves in 3 di- 
mensions. In 2 dimensions the measure is rep- 
resented by a single angle 9 so that in the dis- 
crete representation of a curve, the angle be- 
tween successive segments, A9, varies in the in- 
terval [— 7r,7r]. Unlike the 3d case in which the 
curvature is a positive definite quantity (see Eq. 
[2]), the curvature of a planar curve can be either 
positive or negative (the torsion is zero every- 
where). The direction of the tangent vector can 
be represented as a point on a unit circle and, 
therefore, the measure is given by 

di{s) = d{A9){s) ^ AsrfK(s), in 2d (7) 

We conclude that in summing over the con- 
figurations of 3d curves, one can replace the 
measure in the laboratory frame, D{f{s)}, by 
that in terms of the intrinsic coordinates of 
the deformed line, which is proportional to 
k{s)D{k{s)}. This differs from the correspond- 
ing measure for 2d curve which is proportional 
to oc D{k{s)}. 

3 Curvature distribution of 
worm-like polymers 

In the worm-like chain model (WLC) a polymer 
is modeled as semi-flexible rod whose elastic en- 
ergy is governed only by the curvature 

EwLC = \hl dsK\s) ^ \hY,AsKl (8) 

where h is the bending rigidity. This energy 
is invariant under polar rotations and conse- 
quently the torsion angle Aip is uniformly dis- 
tributed in the range [— 7r,7r]. In 3 dimensions, 
the probability p{K)dn of observing a value of 
the curvature in the interval between k(s) and 
k(s) + dfi; at a point s along the contour of the 
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curve, is proportional to the product of corre- 
sponding Boltzmann factor by the measure ndn 
(see Eq. |6]), 

p{K)dK = /pAse-^'p^^^Vc^K (9) 

where Ip = b/ksT is the persistence length 
(for L » Ip the polymer length does not af- 
fect the statistics of the curvature). Inspec- 
tion of Eq|9] shows that if one rescales the cur- 
vature JlpAs, the probability 
distribution functions of polymers with differ- 
ent bending rigidities can be represented by a 
single master curve ft:exp[— (1/2)k^]. This scal- 
ing may lead to the conclusion that changing 
the persistence length is equivalent to chang- 
ing the discretization. Indeed, calculation of 
properties such as the mean curvature of a 
curve depends on the discretization length As 
and diverges in the continuum limit, (k) = 
\imi\s^o(lpAs)~^^'^ = oo. This, however, is not 
true in general and calculation of physical ob- 
servables such as the tangent-tangent correla- 
tion function, yields well-behaved results: 

(t(0) ■ i{s)) = (to ■ L) = (cos(A^))" (10) 

From Eqsj9]and|5]one can easily conclude that: 
(cos(A6')) = 1 — As/lp ~ exp[— As//p]. Insert- 
ing the above expression into Eq. [TO] one gets 
the WLC relation: 

(t(0)-t») = e-^/'^ (11) 

This correlation function and all properties that 
can be derived from it, such as the mean square 
end to end distance, depend only on the per- 
sistence length and are independent of the dis- 
cretization. 

4 Simulations and analyti- 
cal results 

In order to study the curvature distribution 
function of a worm-like polymer numerically, 
we performed Monte-Carlo simulations of a dis- 
crete polymer made of N equal segments, in 3d. 
To generate the conformations of a polymer we 
used the following moves (see ref. |1]): a point 



between two neighboring segments is picked 
at random and a randomly oriented axis that 
passes through this point is defined. Then, the 
entire part of the polymer that lies on one side 
of the chosen point, is rotated about the axis by 
an angle ip which is uniformly distributed in the 
range [— 7r,7r] (see figure [2]). Since in the WLC 
model the total energy (in dimensionless units, 
ksT = 1) of a discrete polymer is: 

EwLC = bY.[l-cosAen], (12) 

n 

the change of energy in a move is given by 

AEwLc = b [cos(A6'„)oZd - cos(A6'„)„e«;] (13) 

A move is accepted or rejected according to 
Metropolis rule and statistics is collected fol- 
lowing equilibration determined by convergence 
of the total energy (after at least 10^ moves). 
In order to double check our method, we calcu- 
lated the tangent-tangent correlation function 
and the mean square end to end distance as 
well and confirmed that they agree with well- 
known analytical results for worm-like chains 
(not shown). Inspection of Fig. [3] shows that 
there is perfect agreement between our simula- 
tion results and the analytical expression, Eq. 
M The two main features of the distribution 
are the existence of a peak at a finite value of 
the curvature and the fact that the probabil- 
ity to observe zero curvature vanishes identi- 
cally. At first sight, the prediction that any local 
measurement on a worm-like chain will yield a 
non-vanishing curvature, contradicts our intu- 
ition about stiff polymers. Note, however, that 
experiments (e.g., by AFM) that monitor the 
configurations of semi-fiexible polymers such as 
dsDNA, do not directly measure the curvature 
but rather the local bending angle. Even though 
{k) = 1/ ^IpAs does not vanish, the average 
local angle (measured on length scales smaller 
than the persistence length) {66) = ^ As/lp is 
always small. 

The above discussion applies to worm-like poly- 
mers that do not have intrinsic curvature (their 
curvature is generated by thermal fluctuations 
only). Although the WLC model has been ap- 
plied [5] to interpret mechanical experiments 
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on individual double stranded DNA (dsDNA) 
molecules [6], it has been suggested[7] that ds- 
DNA has sequence-dependent intrinsic curva- 
ture. In this case, the experimentally observed 
local curvature of dsDNA contains both the 
contribution of intrinsic curvature and that of 
thermal fluctuations [8]. The generalization to 
the case of spontaneous curvature is straightfor- 
ward: even though the expression for the bend- 
ing energy Eq. [8} is replaced by 
If 1 

(14) 

the measure is not affected and depends only 
on the dimensionality of the problem. Using 
the discrete form of curvature Eq. |5] the distri- 
bution function of curvature for a given intrinsic 
curvature Kq is: 



p{K\Ko)dK = A{KQ)Ke~ 



where 



A{ko) 



J dkKexp{—-l{K — kq 



-1/2 



(15) 



(16) 



is a normalization factor and / is the persistence 
length. An experiment measuring the curva- 
ture distribution of an ensemble of conforma- 
tions of a polymer with a specific primary se- 
quence {ko(s)}, would measure the distribution 
Eq. [15} from which the intrinsic curvature can 
be calculated as a fit parameter. If one moni- 
tors the curvature distribution of an ensemble of 
random copolymers whose sequences are char- 
acterized by a distribution of intrinsic curvature 
p{ko), the corresponding curvature distribution 
function will be: 



dKop{K\Ko)p{KQ) 



(17) 



In order to check this result we performed 
Monte Carlo simulations with different distri- 
butions of ko . As can be observed in Figji], the 
simulations are in perfect agreement with the 
analytical results for the curvature distribution 



experiments. To the best of our knowledge, all 
experiments that probed the local curvature 
of DNA (by electron microscopy [9] and by 
SFM [ini [H]) were done on dsDNA molecules 
deposited on a surface. It has been argued 
that in this case one has to distinguish between 
strong adsorption in which case dsDNA attains 
a purely 2d conformation, and weak adsorption 
in which the 3d character of the adsorbed 
macromolecule is maintained, at least in part 
(see ref. [H]). While in the former case the 
measure of Eq. [7] applies and the most probable 
curvature vanishes [12j, in the latter case one 
should use the 3d measure of Eq. |6] for which 
the most probable curvature is finite and the 
probability to observe zero curvature vanishes. 
New experiments that could probe the 3d 
conformations of DNA and other worm-like 
polymers are clearly necessary to test our 
predictions. 
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calculated using Eq. 17 



5 Discussion 



Finally, we would like to comment on the ap- 
plicability of our results to present and future 
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Fig. 1 3 successive segments are shown by 
black arrows. is the angle be- 
tween the (n — l)'th segment (the ra- 
dius of the unit sphere) and the n'th 
segment. A(y9 is the polar angle de- 
fined by the projection of t„ (dashed 
arrow) on the plane defined by the 
normal ri„_i and the binormal 
(red arrows) 

Fig. 2 Sketch of a Monte Carlo move. The 
right hand part of the curve rotates 
with random angle with respect to 
a random axis (dashed red line) pass- 
ing through a random point m 

Fig. 3 Distribution function of scaled cur- 
vature k for a polymer of = 10^ 
segments of length As = 1 each. 
The empty symbols are the simu- 
lation results for Ip = 25 (black 
squares), Ip = 50 (red circles) and 
Ip = 100 (green triangles). The rigid 
line is the analytical scaled distribu- 
tion function Kexp[— (1/2)^^] where 
k = K^jlpKs. There are no fit pa- 
rameters. 

Fig. 4 In all figures the blue squares and 
green diamonds are the histograms 
of curvature and of intrinsic curva- 
ture, respectively. The red line is the 
analytical distribution function com- 



puted using Eq. 17 
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Figure 1: Caption of Fig. 1. 




Figure 3: Caption of Fig. 3. 
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Figure 4: Caption of Fig. 4. 
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